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Entangling two quantum bits by letting them interact is the crucial requirements for building a 
quantum processor. For qubits based on the spin of the electron, these two-qubit gates are typically 
performed by exchange interaction of the electrons captured in two nearby quantum dots. Since the 
exchange interaction relies on tunneling of the electrons, the range of interaction for conventional 
approaches is severely limited as the tunneling amplitude decays exponentially with the length of the 
tunneling barrier. Here, we present a novel approach to couple two spin qubits via a superconducting 
coupler. In essence, the superconducting coupler provides a tunneling barrier for the electrons which 
can be tuned with exquisite precision. We show that as a result exchange couplings over a distance 
of several microns become realistic, thus enabling flexible designs of multi-qubit systems. 
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I. INTRODUCTION 

Semiconductor based electron spin qubits have made 
significant progress towards scalability. Single qubit gate 
fidelities demonstrated in some devices meet the require¬ 
ments for quantum error correctionii^, with other ap¬ 
proaches not being far behind^^—. Two qubit gates have 
also been realized^ii, and adequate fidelities seem within 
reach^. However, all these gates act over a very limited 
range of typically less than 1 /im, which severely con¬ 
strains scalability as the resulting small qubit spacing 
leaves little room for wiring and local control electron¬ 
ics. For example, a surface code architecture^, which is 
the currently most promising main stream approach to 
error correction, requires a two-dimensional (2D) lattice 
of qubits with nearest neighbor coupling. One possible 
solution is to use charge coupling. Simulations indicate 
that the coupling range can be extended to at least 10 /rm 
with floating electrostatic couplers^ while maintaining 
a strength comparable to that demonstrated in exper¬ 
iments with immediately adjacent qubits^. With this 
coupling strength, the currently achievable level of charge 
noisei^ still leads to coherence times that are two orders 
of magnitude too short to reach the required gate fideli¬ 
ties. Another possibility is to transfer electrons between 
qubits, e.g., using surface acoustic wave o^^d^ or electro¬ 
static gatesi^. First evidence indicates that the spin 
projection can be preserved during such a transfe r^^d*^ . 
While one may hope that one can also achieve spin co¬ 
herent transfer, this remains to be shown experimentally. 
Furthermore, these approaches entail a rather cumber¬ 
some complexity of the device and its operation. 

A very attractive remedy would be to directly extend 
the range of tunnel or exchange coupling between local¬ 
ized electron spins representing a qubit. In principle, this 
could be achieved with a very long and shallow tunnel 
barrier, but in practice disorder would lead to localiza¬ 
tion on a scale of a few 100 nm. Here, we theoretically 
analyze the possibility to use a superconductor that is 
tunnel-coupled to the qubit electrons (Fig. [IJ to medi- 
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FIG. 1. (a) Sketch of the setup analyzed in the paper. Two 
semiconducting quantum dots (white ellipses) are tunnel- 
coupled with strength F to a two-dimensional superconduct¬ 
ing film (dark gray). The levels of the dot are tunable by the 
gate voltages Vi, 2 - (b) Steps in the virtual process leading 
to the exchange interaction. The initial conditions with one 
electron in each quantum dot is shown to the left. Given the 
fact that the electrons are in a singlet state, the following vir¬ 
tual process which is fourth-order in Ht may take place: the 
electron in the left dot tunnels into the superconductor where 
it propagates as a quasiparticle above the gap. The electron 
enters the second dot which becomes doubly occupied. The 
two remaining tunneling processes reverse the path and take 
the electron back to the initial state. Note that that the re¬ 
verse path could be also taking by the other electron which 
leads to a spin flip and explains the factor of two in Eq. 


ate such coupling over extended distances without being 
limited by localization. Qualitatively, the main idea is to 
use extended quasiparticle states for coupling while rely¬ 
ing on the gap to freeze out all the low-energy degrees of 
freedom in the coupler, thus suppressing decoherence. 

A key question is what coupling range and strength can 
be achieved with this approach. To address this question 
for a simple model system, we compute the exchange 
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coupling between two distant electrons (e.g., localized in 
semiconductor quantum dots) that are tunnel coupled to 
a 2D superconducting ground plane and can be detuned 
electrostatically with respect to the latter. We derive 
an expression relating the mediated coupling strength to 
that achievable with direct coupling via the Green’s func¬ 
tion of the superconductor, considering both the ballistic 
and disordered case. Using realistic estimates of the rel¬ 
evant parameters, we find that a micron-scale coupling 
with a practically useful strength of 10 to 100 MHz is 
achievable. 

An important qualitative result is that the decay 
length of the coupling is determined by the detuning be¬ 
tween the quantum dot levels and the upper edge of the 
superconducting gap. This detuning can be controlled 
with high precision via gate voltages. Therefore, the de¬ 
cay length can be substantially larger than the supercon¬ 
ducting coherence length. This result is in contrast to a 
recent proposal!^ considering crossed Andreev Reflection 
as coupling mechanism in a similar setting and finding a 
decay on the scale of the coherence length. 

The remainder of the paper is organized as follows. 
In Sec. m we discuss the general setup and present its 
Hamiltonian. Section Hill discusses the results for the case 
of a clean superconductor. These results are contrasted 
in Sec. m with the results for the disordered case. In 
Sec. El we present realistic experimental parameters that 
allow an exchange coupling over a distance of several mi¬ 
crons. Some technical details about the disorder average 
are moved to the Appendix. 


II. SETUP 


of the system. Note that, if needed, the application of 
a (weak) magnetic field can be included perturbatively 
in the end. We seek a situation where the states |11) at 
energy £i_i = ei -I- £2 and |02) at energy £ 0,2 = 2e2 -f U 2 
are almost degenerate with Se = £ 0,2 — £ 1,1 > 0 much 
smaller than the typical energy spacing in the dots. In 
this situation, we only have to take the states |1,1) and 
|0, 2) into account. The near degeneracy can be achieved 
by setting ei = €2 + U 2 — Se. 

We model the thin film of superconductor by the con¬ 
ventional BCS Hamiltonian Hbcs = 

The spectrum of the superconductor is given by Ek = 
{^1 -\- with A > 0 the energy gap of the supercon¬ 

ductor and ^k = h^k^/2m — here, m is the electron 
mass and p, the chemical potential of the superconduc¬ 
tor. The fermionic Bogoliubov operators Pka are related 
to conventional electronic degrees of freedom via the uni¬ 
tary transformation 

Ckt = Ukl3kt + Vkl3^_k:l, cikl = -'^klSkt + Ukl3^_k:l 


with the parameters Uk,Vk > 0 determined by = 
1 — vl = ^(1 -I- ^k/Ek). In the following, we denote 
with |0) the ground state of the superconductor and cor¬ 
respondingly \ka) = /3fccr|0) denote the (single-particle) 
excitations. 

Coupling between the superconductor and the dot is 
provided by the tunneling Hamiltonian 
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-fH.C., (3) 


We want to study the exchange coupling strength in a 
setup where two semiconducting spin qubits are coupled 
via a thin superconducting film, see Fig. [IKa). The total 
Hamiltonian H — Hd+Hbcs+Ht consists of three parts 
which we will discuss individually in the following. The 
dots can be modeled by the Hamiltonian Hu = Hi + H 2 
with 

~ ^jkij + -^Ujrijfjij — 1 ). ( 1 ) 

It involves the operator rij = Ya counting the 
number of electrons in dot j where nj„ = and 

dja are fermionic operators. The two dots are tunnel- 
coupled via a Hamiltonian to a thin 2D supercon¬ 
ducting film of size L x L. We measure the energies 
relative to the chemical potential of the superconduc¬ 
tor and assume the tunability of the level position e^ 
by nearby gates V). The parameter Uj > 0 describes 
the Coulomb interaction due to the repulsion of multiple 
electrons on a single dot. We note that the eigenenergies 
Erii.na = £ 1*^1 + \Uini{ni - 1) + e 2 n 2 + \U 2 n 2 {n 2 - 1) of 
the dot Hamiltonian Hd only depend on the occupation 
ni,n 2 and not on the spin state of the electrons. This is 
due to the absence of a magnetic field in our description 


where we have taken into account that the two dots are 
separated by a distance R (along the a;-axis). For sim¬ 
plicity, we have assumed t to be the same in the two dots. 
In the following, it will be useful to parameterize t by the 
tunneling rate F/Ii in the normal state with F = 27rt^po; 
here, po = denotes the density of states of the 

normal state (per spin). 

As we are considering the exchange effect mediated by 
the superconductor, an important parameter will be the 
energy difference between the initial state £ip = 2e2 -I- 
E 2 — Se and the intermediate state £ 0,1 + Ek = €2 + Ek 
with the electron in the superconducting wire. For the 
latter to be an excited state, we demand that |£iq — 
£ 0,1 1 < A. In particular, we are interested in a situation 
where £iq — £ 0,1 = €2 + U 2 — Se is smaller but not much 
smaller than A (i.e., the level of dot 1 is tuned close to 
the gap edge). We parameterize this by the energy offset 
M — A — {e 2 + U 2 — Se) > 0 between the initial state with 
one electron in each dot and the intermediate state where 
the electron of dot 1 is transferred to the superconductor. 
Note that M can be tuned independently of Se by £ 2 . 

In this situation, the dominant contribution for the ex¬ 
change comes from a virtual process where we start from 
|1,1) go over to a state |0,1) plus a low energy excitation 
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in the superconductor, then we reach |0,2) before we re¬ 
trace the steps, see Fig.jTJb). In order that this exchange 
interaction can lead to a reduction of the ground state en¬ 
ergy, the spins of the electrons in the initial |1,1) state 
have to be in a singlet as otherwise the |0, 2) state is for¬ 
bidden by Pauli exclusion. The interaction thus assumes 
the form H^x = ■ <^2 with J > 0 the energy dif¬ 

ference between the singlet (which is lowered in energy) 
and the triplet state (which is unchanged). 

Assuming that the spins are in the singlet |l,ls) = 
, we compute the lowering of the 
ground state energy in fourth-order perturbation theory 
in the tunneling Hamiltonian Ht (the triplet energy does 
not change as discussed above). We obtain J = aP^/Jer 
with the dimensionless coupling constant 


r2 


E 

kcr 


(0,2; OlilrlO, 1; ka){0, 1; ka\HT\l, Is; 0) 


£i,i ~ £o,i “ Ek 


1 

27r2p2 
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( 4 ) 


where we have introduced the Green’s function of the 
superconductor 


g{r;E) = 


—i / c?t(0|c^(r)e 
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that we will compute in the following. Note that because 
we are interested in values E = A — M < A, we do 
not need to distinguish between advanced and retarded 
Green’s functions. The exchange interaction is thus given 
by the product of the bare result Jq = T^/fc which would 
be achievable in the case the dots where in direct contact 
and a distance dependent renormalization factor a < 1 
describing the reduction due to the finite spatial separa¬ 
tion. Note that in Eqs. O and we have assumed 
that the superconductor is in the ground state without 
and quasiparticle excitations present which requires that 
the electron temperature is much smaller than the super¬ 
conducting gap A. 


III. CLEAN SC 


In the case of a clean superconductor, it is straight¬ 
forward to evaluate ®. Going over to polar coordinates 
and assuming M <C A <C ^ and kpr 1 yields the 
semiclassical expression 


g{r;E) 


Po dcj) Qi(kF+ik/tiVF)rcos<j> 

^ E-A-ej^A 

^ikpr-iTT I A+i^kT/hvF 

V^27rfc^ y_oo M-|-^^/2A 

/ A X 1/2 

Po cosikpr + 37r/4)e"’'/^^, (6) 


with the effective coherence length ^ = hvp / VSAM that 
is a factor {A/MyP ^ 1 longer than the bare coherence 
length ^0 = Euf/eA. In Eq. (jH), we have taken into 
account that for the relevant part of the integral, we have 
ft! hvpik — kp) ^ A such that u\ ~ \ and Ek ~ 
^ + Cfe/2^- 

Having evaluated the Green’s function in the semiclas¬ 
sical limit, we are in the position to evaluate the dimen¬ 
sionless coupling constant a = 2 cos^ {kpR+2 >p/4) with 


ao 


A 

ApMkpR 
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( 7 ) 


The cos^-dependence of a originates in our model from 
tunnel at point (i.e., momentum-independent tunneling 
amplitudes). In a realistic situation, the diameter d of the 
dot is large such that kpd 1. In this case, the result 
will be modified. In particular, the tunneling amplitude 
will depend on the momentum mismatch between the dot 
and the superconductor. A careful consideration of these 
effects turns out to be rather subtle and beyond the scope 
of the present work, see, e.g.. Ref. [l^. When comparing 
the results for a clean to a disordered superconductor in 
Sec. El we use ao to ease comparison between the clean 
and the disordered case. This corresponds to replacing 
cos^ by its typical value 1/2. In this way, we avoid the 
dependence of the results on microscopic details that are 
relevant only in the clean case. 


IV. DISORDERED SC 


In order to treat the case of a disordered supercon¬ 
ductor it is useful to go over from the Green’s function 
g(r\ E) to the Gorkov Green’s function 


G{r-E) 


d^k (E + efc)e‘^" 

(27r)2 E^-El 


( 8 ) 


as the latter has better analytical properties allowing 
to set up perturbation theory in terms of Eeynman 
diagrams^ii In the limit M <C A <C that we are in¬ 
terested in, we have 


E + ^k _ ul 
E^-El^ E-Ek 


(9) 


such that the two Green’s functions G and g can be used 
interchangeably. It is a well-known resulti^ that under an 
impurity average the Gorkov Green’s function is simply 
given by G{r;E) = G{r; E)e~'^P^ with i the mean free 
path in the disordered system. 

In order to obtain the exchange coupling through a 
disordered superconductor, we should average a over dis¬ 
order. It is important to note that the impurity average 
(denoted by the overline) cannot simply be performed 
separately on the two Green’s functions constituting a. 
The reason is that this neglects interference effects which 
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are relevant for a disordered sample with long phase- 
coherence length. In fact, impurity scattering that in¬ 
volves both Green’s functions even becomes dominant 
at large distances due to the emergence of a new length 
scale, which in the diagrammatic language is subsumed 
by the ladder diagrams forming the diffuson. 

Following ideas of Refs. [T^ and [l^, we calculate the 
diffuson approximation of the product of Green’s func¬ 
tion entering a in the Appendix. We obtain the result 
{E < A) 


J^G{k;E)Gik-q;E)^ 

TTpo 

2 {A^ - E^)[{A^ - E^y/^ + HDq^/2] 


( 10 ) 


with the diffusion constant D = here, G{k;E) 

denotes the Fourier transform of G(r; E). The expression 
m is valid for weak-disorder with kpi ^ 1 and for 
q kp- Using this expression, we can obtain the result 
for the disorder-averaged exchange coupling {q, (j) are the 
polar coordinates of q) 


A f (Pq (,igRcos^ 

SttpoM J (27r)2 y/2AM + hDq^/2 

AKo{R/^d) ^ , 

2TTMkFe ^ 2{2nyGMkF^R^^‘^^ 


( 11 ) 



FIG. 2. Plot of the dimensionless coupling constant for clean 
aluminum (oq) and disordered aluminum (a) as a function of 
the distance R of the spin qubits that are coupled. The plot 
shows both the results for M/A — 5 x 10“® (solid lines) and 
the results for M/A = 5 x 10“"* (dashed lines). 


gap 7 (the so-called Dynes parameter) should be much 
smaller than M. Recent experiments have shown that 7 
can be as small as 10“®A^ Given this input, we take a 
conservative choice of M = 1 peV which corresponds to 
M/A = 5 X 10“^. As a result, we obtain the clean effec¬ 
tive coherence length ^ = 36 pm and the dimensionless 
coupling constant ao of Eq. ©, see Fig.H 


with 

When comparing the ballistic result © to the diffu¬ 
sive result (El) there are two major differences: (i) the 
exponential decay is controlled by different length scales 
^ versus ^f>. (ii) the algebraic decay of the former is 
given by R~^ whereas the latter decays more slowly with 
the power 

V. ESTIMATE OF THE COUPLING 
STRENGTH 

We would like to end by discussing realistic length 
scales R over that the superconducting film can be em¬ 
ployed in order to exchange couple two spin qubits. Start¬ 
ing from realistic values Jg — fi(10 — 100) GHz for the di¬ 
rect exchange coupling of two qubits, we aim at achiev¬ 
ing a dimensionless coupling constant a ~ 10 “^ in or¬ 
der to end up with useful exchange coupling constants 
J ~ /i(10 — 100) MHz. For the superconductor we pro¬ 
pose aluminum with a gap parameter A/ks = 2 . 2 K 
which corresponds to a coherence length ^0 = 2.3 pm 
for a clean sample. Aluminum has a Fermi velocity 
vf = 2.0 X 10® m/s that implies a Fermi wave-vector 
kp = ^7 nm“^. The most crucial parameter which makes 
it possible to increase a is the detuning M. Some of the 
requirements bounding M from below is that the detun¬ 
ing should be held stable over the time of exchange in¬ 
teraction and that the smearing of the superconducting 


ao 


0.94 nm 
R 


— i?/36 
^ 5 


( 12 ) 


in a clean system as it is for example obtained by epitax¬ 
ial growthi^ For a distance of i? = 1 ^m this evaluates 
to ao = 9.1 X 10“"^. We note that the value of ao is 
completely dominated by the prefactor. Thus, we expect 
that the analysis becomes more favorable for the case of 
disordered aluminum as is obtained for example by sput¬ 
tering. 

Extrapolating from Refs. [23l - I^ . a realistic value of the 
mean free path in aluminum is ^ = 100 nm which trans¬ 
lates to a diffusive coherence-length = 1.3/rm. This 
yields 


S= (13) 

and evaluates for i? = 1 pm to the more favorable result 
a = 1.3 X 10-2. 

Figure [5] shows that for distances R smaller than a 
characteristic distance R* the coupling via a diffusive 
superconductor is larger than the one through a clean 
system^S. For R> R* this reverses. For the typical case 
£ <C the characteristic distance R* is approximately 
given by 


R* 


2 


ln(7r^/2€). 


(14) 


For M/A = 5x 10 we find numerically that R* = 
5.9 pm and that a is larger than 10“® up to i? = 3.6/rm. 
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Provided accurate control over the gate voltages, the 
range of the exchange coupling can be extended even fur¬ 
ther. For example, lowering M by a factor of 10, i.e., for 
M = 0.1 ^eV we obtain the dashed lines in Fig. O We 
see that for distances smaller than R* — 12 fim the dif¬ 
fusive material leads to a stronger exchange interaction. 
Furthermore, a is larger than 10“^ for distances up to 
i? = 11 /im. 


VI. CONCLUSION 

We have derived the strength of the exchange inter¬ 
action between two spin qubits which are connected via 
a 2D him of superconducting material acting as a cou¬ 
pler. We have shown that the bare exchange interaction 
for two neighboring dots is reduced by a dimensionless 
coupling constant a incorporating all effects of the hnite 
distance R. We have presented results for the case of 
both a clean and a disordered superconductor. We have 
shown that for distances R smaller than a characteris¬ 
tic distance i?*, a diffusive superconductor outperforms 
the clean one due to the prefactor in a for the former 
being a factor (RjlY!'^ larger than for the latter. We 
have shown that a diffusive superconductor with a mod¬ 
erate mean free path of ^ = 100 nm enables to have use¬ 
ful exchange coupling strengths over a distance of more 
thant 10 /rm. In practice, it would be convenient to use 
superconducting wires rather than an extended film to 
mediate the coupling. In this case the further confine¬ 
ment would eliminate the decay of the prefactor in 

Eq. (fTTl) . leaving only the exponential decay and allow¬ 
ing substantially stronger coupling. This suggests that 
superconducting wires might be suitable for mediating 
an exchange coupling between quantum dots with a sep¬ 
aration large enough to implement a 2D lattice. An ad¬ 
ditional advantage is that the coupling strength can be 
varied easily and rapidly by changing the electrostatic 
potential of the superconductor. 


ACKNOWLEDGMENTS 

FH acknowledges financial support from the Alexander 
von Humboldt-Stiftung. 


Appendix A: Calculation of the Diffuson 

In this Appendix, we apply the ideas of Refs. [I^ and 
[H to the specific case of disordered 2D superconducting 
film. The calculation of the impurity average has to be 
performed in Nambu space, because the Green’s function 
has anomalous components due to the superconducting 
condensate. All the Green’s function can be combined in 
the matrix Green’s function 

nil, E + ^kT 3 + ATi 

Gik] E) = - E^-El - 


acting via the Pauli matrices Ti on the Nambu space. 
Note that the connection to the Gorkov Green’s func¬ 
tion introduced in the main text is given by G(fc; E) = 
G{k,E)ii. We assume throughout this section that 
0 < A < A such that the excitations are virtual and 
similar to the Matsubara formalism we do not have to 
worry about retarded and advanced Green’s function. If 
required, the results for A > A can simply be obtained 
by analytical continuation. 

We are interested in performing an impurity average 
in the potential V{r) = ~ u)'r 3 where 

denotes the position of the i-th impurity. The connec¬ 
tion with the mean-free path is given by the Born result 
Hyp/i = 2'KniV ^pq where rii is the density of impurities. 
The impurity averaged Green’s function is given by (see 

Ref. 113) 


G{k-E) 


E + ikT3 + Ari 

E^-e^- 


where E = pE, A = ryA, and 

=1 , 

2£(A2-£;2)i/2' 


(A2) 


(A3) 


In order to calculate the disorder average of a product 
of Green’s functions, it is useful to introduce the four 
matrix diffusons (tq is the identity matrix) 


DM = J ^^3G(fe; A)T,G(fc - q; E^ (A4) 

In terms of these, the combination of Green’s functions 
in m is given by 

J^G{k; E)G{k -q;E) = ^[Doiq) + DsiqM 

. 

The dominant contribution to in the weak disorder 
limit originates from ladder diagrams which keep the rel¬ 
ative momentum q conserved. The zeroth order term is 
given by 

Dl°\q) = nv^ J^r3G{k-E)nG{k-q-EM (A6) 

As we are interested in small relative momenta q, we can 
expand and obtain D^^\q) = ai + biq^ with 

2 f ,1k A^ - AEti 

ao = nv poj d^k T 3 G{k- E) rg = 

(A7) 


The expansion in q is obtained by the replacement 
= ^k — kvpq cos 4> with (f> the angle between k and 
q. The first order in q vanishes due to the integration 
over (j). The first non-vanishing contribution reads 

fi^Vpai 

" ~ 8(A2-A2)' 


(AS) 
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As has a term proportional to ri, we need to calcu¬ 
late additionally in order to obtain a closed set of 
equations. We obtain the expression 


hvp AA — E'^ti 


(A9) 


the linear system of equations, we obtain the result 


A(g) = 


I?, 


( 0 ) 




i(o) 




(A12) 

Inserting the expressions for yields the final result 


For completeness, we give also the other expressions 


hvp T2 
'2f (A2 -.^2)1/2’ 


03 = 0. (AlO) 


For further convenience, we denote by Df'j the term in 

Dproportional to tj such that by 

definition. 

Summing the ladder diagrams is equivalent the Dyson 
type equations^! 


A(g) = Df\q)+Y,D^^]{q)D,{q)- (AH) 

3 


The system is closed in the subspace i,j G {0,1}. Solving 


Doiq) = 


hvp 


A2 - AFln 


2£ (A2-F;2)[(A2-F;2)i/2 + 7i£)g2/2] 

(A13) 

valid for small q. For reference, we also give the result 


Di{q) = 


hvf 


AE - E'^ti 


21 (A2-F;2)[(A2-F;2)i/2 + ;i£)g2/2]- 

(A14) 

The diffuson only couples to itself, thus the solution 
of (|A11I) is simply 


n r 1 = ^ _72_ 

" 1-D® 2£ {A'^-E'^Y/^ + hDq^/2 

(A15) 

The Z ?3 diffuson vanishes such that only Dq enters the 
expression jAH). From (IA5I1 , we obtain the result quoted 
in the main text. 
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